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To study the critical Casimir force between chemically structured boundaries immersed in a binary 
mixture at its demixing transition, we consider a strip of Ising spins subject to alternating fixed spin 
boundary conditions. The system exhibits a boundary induced phase transition as function of the 
relative amount of up and down boundary spins. This transition is associated with a sign change of 
the asymptotic force and a diverging correlation length that sets the scale for the crossover between 
different universal force amplitudes. Using conformal field theory and a mapping to Majorana 
fermions, we obtain the universal scaling function of this crossover, and the force at short distances. 

PACS numbers: 11.25.Hf, 05.40.-a, 68.35.Rh 


Fluctuation-induced forces are generic to all situations 
where fluctuations of a medium or field are confined by 
boundaries. Examples include QED Casimir forces mm, 
van der Walls forces [3], and thermal Casimir forces in 
soft matter which are most pronounced near a critical 
point where correlation lengths are large W- The inter¬ 
action is then referred to as critical Casimir force (CCF). 
Analogies and differences between these variants of the 
common underlying effect have been reviewed in Ref. [B]. 

Experimentally, CCFs can be observed indirectly in 
wetting films of critical fluids [7], as has been demon¬ 
strated close to the superfluid transition of 4 He |S] and 
binary liquid mixtures [5j. More recently, the CCF 
between colloidal particles and a planar substrate has 
been measured directly in a critical binary liquid mix¬ 
ture mm- Motivated by the possibility that the lipid 
mixtures composing biological membranes are poised at 
criticality Buna, it has been also proposed that inhomo¬ 
geneities on such membranes are subject to a CCF [Tl] 
which provides an example of a 2D realisation. 

The amplitude of the CCF is in general a universal scal¬ 
ing function that is determined by the universality classes 
of the fluctuating medium m- It depends on macro¬ 
scopic properties such as the surface distance, shape and 
boundary conditions of the surfaces but is independent of 
microscopic details of the system [5j. Controlling the sign 
of fluctuation forces (attractive or repulsive) is important 
to a myriad of applications in design and manipulation 
of micron scale devices. While for QED Casimir forces 
a generalized Earnshaw’s theorem rules out the possibil¬ 
ity of stable levitation (and consequently force reversals) 
in most cases m, the sign of the CCF depends on the 
boundary conditions at the confinement. For classical bi¬ 
nary mixtures, surfaces have a preference for one of the 
two components, corresponding to fixed spin boundary 
conditions (+ or —) in the corresponding Ising universal¬ 
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FIG. 1: Ising strip of width L with alternating fixed spin 
boundary conditions on one side, with a typical spin con- 
figuation indicated by the shading. 

ity class. Depending on whether the conditions are like 

(++ or-) or unlike (H— or —h) on two surfaces, the 

CCF between them is attractive or repulsive. So-called 
ordinary or free spin boundary conditions are difficult 
to realize experimentally but can emerge due to renor¬ 
malization of inhomogeneous conditions as we shall show 
below [17] , 

Motivated by their potential relevance to nano-scale 
devices, fluctuation forces in the presence of geomet¬ 
rically or chemically structured surfaces have been at 
the focus recently. Sign changes of CCFs due to wedge 
like surface structures have been reported very recently 
[Tgj . Competing boundary conditions can give rise to in¬ 
teresting crossover effects with respect to strength and 
even sign of the forces. Here we consider such a situ¬ 
ation for the Ising universality class in 2D. At critical¬ 
ity, this system can be described by conformal field the¬ 
ory (CFT) and CCFs are related to the central 

charge of the CFT [2Tfj2fT| . and scaling dimensions of 
boundary operators [231 . 

In this Letter, we show that boundary conditions which 
alternate periodically between two spin states (see Fig. [l]) 
give rise to a novel phase transition. Associated with that 
is a diverging correlation length that sets the scale for a 
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sign change of the CCF on one side of the transition. We 
obtain the critical exponents and exact expressions for 
the universal scaling function of the force in the critical 
region. Consider the Ising model on an infinitely long 
strip of width L, and assume that the system is at its 
critical temperature T c so that it is conformally invari¬ 
ant. For homogenous, fixed spin boundary conditions 
71 , 72 = ± on the two boundaries, the critical Casimir 
energy per unit strip length, T , is determined by CFT. 
Since L is the only finite length scale, the energy obeys 
a simple power law. The amplitude is determined only 
by the central charge c = 1/2 of the Ising model and the 
scaling dimension /i 7l72 of the so-called boundary condi¬ 
tion changing (BCC) operator from 71 to 72 (see below 
for details on the BCC operator) [23] , 

■ F= -7F (ig “ ^7172^ J (!) 

where we measure here and the following energies in units 
of fcflT c . For like boundary conditions 7i = 72 = + 
or — one has /i++ = h _= 0 and hence an attrac¬ 

tive force F = —dT/dL. For unlike boundary conditions 

71 ^ 72 one gets h -1 _= h _j. = 1/2 and hence a repulsive 

force. There is one more conformally invariant boundary 
condition that corresponds to free (f) spins or ordinary 
boundary conditions. When combined with fixed bound¬ 
ary conditions, the corresponding BCC operator has the 
scaling dimension hf + = hf _ = 1/16 which implies a 
repulsive interaction in Eq. |l]). In the following we con¬ 
sider a strip with homogeneous + spins on one boundary 
and alternating regions of — and + spins of length a and 
b, respectively, on the other boundary, see Fig. [T] 

If the temperature is slightly different from T c , the 
system is in the critical region, where the free energy 
density can be decomposed into non-singular (F ns ) and 
singular (F s ) contributions, 

F(t, L, t) = F ns {t , L, r) + F s (t, L, r) (2) 

that depend on the reduced temperature t = T/T c — 1, 
the width L, and a scaling variable r = a/b — 1 that is 
specific to the alternating boundary conditions in Fig. |Tj 
While the non-singular part is an analytic function of t 
and r, the singular part is not. For homogeneous bound¬ 
ary conditions, t is the only relevant scaling variable, and 
in the critical region the singular part of the free energy 
density is given by a universal scaling function d that 
depends only on L/£ [5l ITTTj where £(i —> 0 ± ) = 
is the bulk correlation length with amplitude and ex¬ 
ponent v = 1 for the Ising model. As we shall see below, 
the same renormalization-group (RG) concepts apply to 
a novel, boundary induced critical region that we identify 
for inhomogeneous boundary conditions around a = b. 
To focus on that region, we assume in the following that 
the system is at its bulk critical point, t = 0. For large 
L a,b the singular part of the free energy density can 


be expressed in terms of a universal scaling function of 
the new correlation length £ c (r) = (o + &)|t| _ " 0 , 

7/(0 ,L,t) = ±0[L/&t)]. (3) 

Below we shall determine d and the exponent v c . 

BCC operators have been introduced in CFT to study 
systems with discontinuous boundary conditions 24] . 
When inserted on a boundary, these local operators in¬ 
terpolate between the different boundary conditions on 
either side of the insertion point. They are highest weight 
states of weight h and all such states may be realized by 
an appropriate pair of boundary conditions. For the crit¬ 
ical Ising model, the BCC operator that takes the bound¬ 
ary condition from + spin to — spin corresponds to the 
chiral part of the energy operator e(z, z). This can be un¬ 
derstood easily in the representation of the Ising model in 
terms of a free Majorana fermion field ip(z) out of which 
the energy operator is composed, e(z, z) = iil>(z)ijj(z) 
[25l : The Jordan-Wigner transformation shows that the 
fermion creation and annihilation operators flip locally 
the spin orientation. 

Now the BCC operators permit us to relate the parti¬ 
tion function of the strip with alternating boundary con¬ 
ditions to a correlator for the field ’ip(z) at positions where 
the boundary conditions change. On the upper complex 
plane, one has ( tp(z)ip(z ')) = 1 /{z — z') which yields (af¬ 
ter a conformal map) for the partition function of the 
strip the Pfaffian, 

Z = ZoMwi) • = Z 0 Pi(G) = Z 0 det 1 / 2 (G), 

(4) 

with G = [(ip(wi)'4 , (wj))]i t j=i : ... t 2Ni where we used the 
Wick theorem for fermions, Wj are the positions of the 
2N BCC operators on the upper edge of the strip, and 
Zq is the partition function of the homogenous system 
with a = 0. Due to the symmetry under translations by 
a + b , the matrix G is of block Toeplitz form, G t j = gi~j, 
with 

n ._( g\j(a + b )] g\j{a + b)-a]\ , , 

9j V3[-?'( a + & )+“] 9li(a + b)\ )' 


where g(w) = n/[2L sinh( 7 rui/( 2 L)]. 

The free energy density can be expressed in the ther¬ 
modynamic limit as 


T = - 


7T 1 

48 Z 


1 


iV -400 2IV(a -(- b ) 


log det G. 


( 6 ) 


The Szego-Widom (SW) theorem for block Toeplitz ma¬ 
trices states that the determinant can be expressed 
in terms of the matrix valued Fourier series tp(9) = 
Er=-oo 9ke ike as (26; 


lim —— log det G 
00 2N 6 



dfllogdet v?((?) (7) 
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where det acts now on a 2 x 2 matrix. It turns out that 
this formula can be only applied for the case a < b. The 
reason for that is a subtle difference between the Toeplitz 
matrix G and the corresponding circulant matrix C that 
describes periodic boundary conditions along the strip. 
While for a < b the spectra of G and C become equivalent 
for N —> oo, for a > b there exists a pair of eigenvalues 
of GC _1 that tend to zero exponentially for N —> oo, 
yielding an extra contribution S that is determined by 
the decay of the Fourier integral 


1 f 2n 

J = — y d9e~ lJ0 [^- 1 (0)] 11 ~ e~ jS for j -» oo (8) 


and has to be subtracted from the r.h.s. of Eq. ([7]) for 
a > b. Here [^ _1 (0)] denotes the 11-element of the 
2x2 matrix c^ -1 (0). In the following we apply Eqs. ([T]) 
and (|8j) to compute the critical Casimir force in various 
scaling limits. 

When L <C a, b, the function g(w) defined below 
Eq. jH]) can be replaced by g(w) = (n/L)e~ n ^^ 2L ' > 
which yields the exact determinant 


det <p(0) 


cos 9 — cosli(7r(a — b )/(2 L)) 
cos 9 — cosli(7r(a + b )/(2 L)) 


(9) 


For a < b, the SW theorem then yields ^ log det G = 
— (7ra)/(2L). For a > b, this is also the correct result as 
it follows from subtracting the correction S which follows 
from Eq. ^ and J = e -’ r l°- 6 b'/( 2i ) as 5 = w(a — b)/(2L). 
It follows that the critical Casimir force for L <C a, b is 


7 r 23a — b 1 
48 a + b L2 + '" 


( 10 ) 


It has an analytic amplitude that varies continously with 
a/b. This result is identical to an addition of the am¬ 
plitudes from Eq. ([!]) for unlike and like boundary con¬ 
ditions, weighted by a/(a + b) and b/(a + 6), according 
to their occurrence. Hence, additivity holds at short dis¬ 
tances. This has been observed also for a 3D Ising model 
in the special case of boundaries with alternating stripes 
of equal width [Hi- 

Next, we consider the case L a,b. Using the Abel- 
Plana summation formula, it can be shown that in this 
limit the elements of the matrix 


**) = z 


( *7i(6>) 

1 - 75 ( 0 ) 


72 ( 9 )\ 
*7i(0)y 


( 11 ) 


approach 


7 j(9) = { 7 j{9, t) — i J+1 [tanh(0L/(a + b)) + tanh((0 — 2n)L/(a + 6))]} 


( 12 ) 


with 71 (9,t) = 1 — 9/n and 


72 (9,t) 


1 

7r 


t + 2 

T + 1 


+ e l \r + 2) 2 F l 



1 t + 3 
r+ 2’ r + 2’ 



— e 


-ie 


t + 2 
2t + 3 


2-Zl 


/ 2r + 3 3r + 5 
V 1, T + 2 ’ t + 2 



(13) 


where 2 E 1 is a hypergeometric function. For a < b, the SW theorem then yields the free energy density 


T=- 


7r 1 

48 Z 


1 

47 r(a + b) 



1 + F( 0 ,r) 


tanh 


9L 
a + b 


+ tanh 


(9 — 2n)L 
a + b 


d9. 


(14) 


where we have subtracted a L-independent contribution 
that does not change the force, and defined T( 0 ,r) = 
(2yi + i (72 - 72 ))/(l 7 i| 2 - |T 2 1 2 )- In the evaluation of 
the integral, the correlation length £ c (r) defined above 
Eq. § becomes important. The integrand is exponen¬ 
tially localized around 9 = 0,2n over a small range 
(a + b)/L. Also, it can be shown that T has the scaling 
property lim T ^ 0 T(r 2 /C, r) = r 0 (C) = 1/(1+tt 3 C/32) for 
any constant Hence, in the critical region of small r, 
or £c( t ) a + b, the proper scaling is obtained by setting 
C = (Tr 2 )/(a + 6 ) = L/£ c (up to a numerical coefficient), 
showing that the exponent v c = 2. In the integral, T(f?, r) 
can be replaced by To(C) and one obtains after a simple 
integration the result for the universal scaling function of 
Eq. § when a < b or t < 0, 




- 1 


(15) 


where Li 2 (x) = Y2kLi xk /^ 2 ^ a polylogarithm func¬ 
tion. Outside the critical region L £ c , one has 
$_(£ —> 00 ) = —7t/48 so that the force is fully domi¬ 
nated by the boundary regions with like spins. On the 
contrary, for L -C £ c , and hence r —► 0 _ , the frustration 
between almost equal amounts of fixed + and — spins on 
the boundaries leads to a renormalization to effectively 
free boundary conditions with d-(( —> 0) = 7t/24. For 
a > b, the correction S yields an extra contribution Ad(Q 
determined by 

Atf (C) tan[Atf (C)] = (16) 

so that the scaling function for r > 0 is $+(C) = 
d-(C) + Ai?(£). Since Ad(( —► 0) = 0, the scaling func¬ 
tion is continuous around r = 0. For L £ c , how- 


1 + tt 3 C/32 
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ever, Ad(£ —► oo) = 7 t/2 so that the system asymp¬ 
totically realizes homogenous unlike boundary conditions 
with $+((, —y oo) = 237r/48. 



23 fc 

FIG. 2: Schematic overview of critical Casimir force ampli¬ 
tudes as function of the strip width L and the ratio a/b. For 
L 3> a, b the solid curves represent the diverging correlation 
length £ c . The horizontal dashed line indicates the cut along 
which the force amplitude is plotted in Fig. [3] Along the red 
curve the sign of the force changes whereas the blue curve in¬ 
dicates only a change between two universal (repulsive) limits. 


Our findings can be summarized by the scheme of 
Fig- [?} It shows the different scaling regimes and the cor¬ 
responding asymptotic amplitudes of the Casimir force. 
At short distance L <C a, b the amplitude varies continu¬ 
ously across the critical point at a = b, with a sign change 
at b/a = 23. For L a, b there exist three distinct re¬ 
gions: around a = b appears a region where L <C £ c where 
the force is repulsive and approaches for asymptotic L the 
universal amplitude for fixed-free spin boundary condi¬ 
tions. For a < b, the force changes sign from attractive 
to repulsive when L approaches £ c , corresponding to a 
stable point. For a > 6 , the force is always repulsive but 
the amplitude crosses over from 7 t/ 24 to 237 t/ 48 under 
an increase of L beyond £ c . 


The dependence of the force F on \a — b\ at fixed 
L a, b (see dashed horizontal line in Fig. [2]) is de¬ 
termined by F = —dF/dL = 0(x S )L ~ 2 with a uni¬ 
versal scaling function 0 of the scaling variable x s that 
is defined on both sides of the critical point by x s = 
sign(r)(L /^ c ) 1 / 2 ~ a — b. This function is shown in Fig. 3] 
where we used the results for d±(A/£ c ) of Eqs. (15), (16). 



-5 0 5 10 

j: s =sgn(a-b)(L/f f ) 1,2 =[L/(a+b)] l ' 2 (a/b-l) 


FIG. 3: Universal scaling function 0(x s ) for the critical 
Casimir force as function of the scaling variable x s ~ a — b. 


In the critical region |cc s | <C 1, one has the expansions 


0(x 5 ) = 


- m x s + ■ ■ • for x s < 0 

3/2 

+ § 77 x s + ■■■ for x s > 0 


(17) 


whereas for L outside the critical region, |x s | 1, 

for x s < 0 


0(x s ) = 


7r | 32 log 2 1 
" 48 ' 7r 4 xl 


237T _ 32(7r 2 —log 2) J_ for X > 0 

48 7T 4 x 2 + ■ ' ' IOT Xs > U 


• ( 18 ) 


We see that <d(x s ) is not analytic around x s = 0 and 
hence constitutes the singular part of the free energy 
density, see Eq. ([2|. This resembles the singular na¬ 
ture of scaling functions describing the bulk transition 
at T = T c . 

Our results show the existence of a novel phase transi¬ 
tion for the critical Casimir force in the 2D Ising model 
that is induced by inhomogeneous boundary conditions 
with a varying ratio of up and down spins. We ob¬ 
tained exact expressions for the universal scaling func¬ 
tion of the force. Due to the observed renormalization of 
boundary conditions, in binary mixtures, ordinary (free 
spin) boundary conditions can be realized experimentally 
and “switched” on and off by varying the distance L, or 
an inhomogeneous surface field. The crossover between 
different universal amplitudes leads to a stable equilib¬ 
rium point at L ~ £ c for 1/23 < a/b < 1. The emer¬ 
gence of the novel phase transition at a = b is related 
to the relevance of a surface magnetic field ~ r at a 
surface with free spin boundary conditions. This can 
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be seen from the decay of the spin correlations along a 
single surface, ( a x a x /) ~ \x — with = 1 for 

free boundary conditions [I5][27]. Since the surface field 
contributes an energy ~ r f dxa x , the scaling dimension 
y T = l/z/ c = 1 — 77||/2 which is identical to our find¬ 
ings above. It is interesting to explore these concepts in 
general spatial dimensions for Ising and XY models, and 
tri-critical points which have an even richer spectrum of 
possible boundary conditions. 

We thank M. Kardar for many fruitful discussions. 
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